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Abstract 



We present a quenched lattice calculation of heavy-light meson 
decay constants, using non-relativistic (NRQCD) heavy quarks in the 
mass region of the b quark and heavier, and clover-improved light 
quarks. The NRQCD Hamiltonian and the heavy-light current include 
the corrections at first order in the expansion in the inverse heavy 
quark mass. We study the dependence of the decay constants on the 
heavy meson mass M, for light quarks with the tree level (csw = l)j 
as well as the tadpole improved clover coefficient. We compare decay 
constants from NRQCD with results from clover (csw = 1) heavy 
quarks. 

Having calculated the current renormalisation constant Za in one- 
loop perturbation theory, we demonstrate how the heavy mass de- 
pendence of the pseudoscalar decay constants changes after renormal- 
isation. For the first time, we quote a result for fs from NRQCD 
including the full one-loop matching factors at 0(a/M). 



PACS numbers: 12.38. Gc, 12.39.Hg (HQET) 13.20.He (leptonic and semilep- 
tonic decays of bottom mesons), 14.40.Nd 



1 Introduction 



A calculation of the decay constant of the B meson, is of interest for 
the determination of the unitarity triangle parameterising CP-violation in 
the Standard Model. The element \Vtd\ of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix can be determined from an experimental study of B — B 
mixing, using f B as one of the input parameters |I[ . It can be defined through 
the following matrix element (in Minkowski space): 

ipJ B = (OlA^B), (1) 

where is the heavy-light axial vector current. Similarly, the vector decay 
constant can be defined through the relation^: 

ie„f v M v = (0\V„\B*), (2) 

being the heavy-light vector current. A practical tool for calculations 
involving hadrons containing one heavy quark is Heavy Quark Effective The- 
ory (HQET) ||. It exploits the fact that in heavy-light systems, in the limit 
of infinite heavy quark mass itlq, there is a spin-flavour symmetry between 
heavy quarks. Corrections due to finiteness of the heavy quark mass are in- 
cluded in an expansion in l/m,Q. For the decay constant one expects a heavy 
mass dependence of the following kind: 

fVM = Ao(l + At/rag + A 2 /m 2 Q + ...). (3) 

Here, / denotes the pseudoscalar or vector decay constant and M the mass 
of the corresponding heavy meson. Using lattice QCD, these matrix elements 
can be calculated nonperturbatively from first principles. Lattice calculations 
of heavy-light decay constants have been performed using the relativistic 
'naive' (see e.g. ||), or the clover-improved (e.g. Wilson action for the 
heavy quarks. A reinterpretation of the naive relativistic action in the regime 
arriQ > 1 has been proposed by 0. In recent studies (e.g. || [F], |8|, 0), 
this suggestion has been implemented to various degrees. These simulations 
have either used heavy quarks in the mass region of the charm, relying on 
extrapolations to the b quark, or gone at most up to masses around the b. 
Lattice calculations have also been done in the infinite mass (static) limit (e.g. 

1 With this definition /y = fs in the static limit 
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Refs. 


10 
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see Ref. 
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121). For a recent overview of the status of lattice calculations 



Alternatively, heavy quarks can be simulated using non-relativistic QCD 
(NRQCD) fl4]| , an effective theory where the operators in the action and 
heavy-light currents are expanded in a series in the bare inverse heavy quark 
mass rrvq . With this approach, one can study quarks in the whole region 
between the b quark and the static limit. The first calculation of heavy- 
light decay constants with NRQCD [[15], and a following more extensive 



simulation [[16] used an NRQCD action at 0(l/rriQ ) for the heavy quark 

and a Wilson action for the light quark; however the 0(l/mg' ) ) corrections 
to the current operators were not included. A calculation using quenched 
configurations and Wilson light quarks that includes the current corrections 
is introduced in Ref. [17]. The first report on a study of decay constants where 
the currents are also corrected through 0{l/rriQ > ), using Wilson light quarks 



and dynamical configurations, is published in fL8|| . In the study described 
in this paper, we use an NRQCD action and currents as in [TJ|] but with 
quenched configurations and clover light quarks. 

In Sec. |2|we explain the operators that we use in our NRQCD Hamiltonian 
and heavy-light currents, the simulation parameters and the interpolating 
operators for the mesons. The fitting procedure and results for energies and 
amplitudes are presented in Sec. |3|. The analysis of the results follows in 
Sec. (|. After a brief discussion of the heavy-light meson masses using tree- 
level- improved (csw — 1); an d tadpole-improved light fermions, we turn to 
the decay matrix elements. We compare the decay constants with tree-level- 
improved light quarks and NRQCD heavy quarks with results using a clover 
action with csw — 1 a ls° for the heavy quark. This is followed by a study of 
axial, vector, and spin-averaged matrix elements with tree-level and tadpole- 
improved light clover quarks as a function of the heavy quark mass. We 
give results for the physical ratio fs a / fB d - Using a renormalisation constant 
Za from one-loop perturbation theory, we finally present renormalised axial 
matrix elements, and quote our estimate of fs- 

2 Simulation details 

We choose to work in a Pauli basis where the two-component heavy quark 
spinor Q and the antiquark spinor decouple |[18||. Only the spinor Q is 



3 



used in our simulations. The non-relativistic Lagrangian which describes the 
b quark is expanded through O^/m^g) at tree level: 

C = QHD t + H + 6H)Q, (4) 

where 

The gauge links are tadpole improved: 

U» -> C^/«o, ^ = (l/3Tr£/ Plaq .>, (6) 

so the covariant derivatives act in the following way on a Green function G: 

A^G(x) = [U fM (x)G(x + f l )-Ul(x-fi)G(x-il)}/(2u ); (7) 

A®G(x) = [U^(x)G(x + p)+Ul(x-fi)G(x-fi)]/u -2G(x), (8) 

and 

A (2) =EAf (9) 

We use the tree level coefficient cb = 1 in Eq. Using tadpole-improvement, 
we expect that the perturbative contributions to this coefficient do not be- 
come too large. For q* values ranging between 1/a and ir/a, ay lies between 
~0.25 and ~0.15. The perturbative correction to this operator may thus be 
roughly of order 20%. The heavy quark follows the evolution equation [|l9j 

on the first timeslice, and on the following timeslices 

aH \ n aH^ n 



Gt+1 = ( 1 ~ ~2^J Ul V ~~ "2^ ' (1 " a5H ^ 1 > L ( U ) 



It has been noted |2(J that this evolution equation introduces an error in 
principle of 0{aK^ CD /m ( Q > ) in the amplitude since the operator 5H is not 
applied on the first timeslice. We estimated the actual size of the error 
by comparing it with a different evolution equation and found it to be ~ 
3 — 4% for the bare lattice matrix element plj (using clover light quarks and 
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dynamical configurations). This deviation is of the order of the statistical 
error (see Sec. (|), so it will be ignored here. 

In the calculation of decay constants to the desired order in the 1 / Mq'- 
expansion, one has to also include the corrections to the currents. At tree 
level, these can be obtained by relating the heavy quark field in full QCD, h, 
and the nonrelativistic heavy quark field through an inverse Foldy-Wouthuysen 
transformation on the heavy quark spinor. At 0(l/mg' ) ), one has: 



)-*«■»)(«), (12) 
where 



We write the heavy-light currents in full QCD as 

Aft = <?757aA and = q^^h, (14) 

q being the light quark field. In the following we consider only the time 
component of the axial vector current and the spatial components of the 
vector current. In our simulations the conventions for the 7 matrices are: 

70 = { -1 J ' 7 = { -iff J ' 75 = { 1 J ' (15) 

The current, corrected though 0(l/mg' ) ), takes the form, 

J» = jf ] + J£\ (16) 
where the contributions to the currents are given by: 

4 0) = QuQ, JP = -i-^m^y ■ DQ (17) 

2m Q 

for the axial current and 

4 0) = -iqt^JkQ, 4 1J = ^<?W ■ 3 Q ( 18 ) 
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for the vector current, where we use the notations qu for the upper and ^34 
for the lower two components of the light quark spinor. Other operators of 
the same mass dimension and lattice symmetry mix under renormalisation 
with the Foldy-Wouthuysen operators. For the axial current, we discuss this 
further in subsection |4.3| . A more general list of these operators can be found 
e.g. in Refs. |2f and |j§. 



We also implement the heavy quark in the static approximation, which 
corresponds to the Lagrangian: 

C = Q ] D t Q. (19) 

The static heavy quark propagator follows the evolution equation: 

G t+l - U\G t = 5 X , . (20) 

Our light quark propagators were generated by the UKQCD Collaboration. 
These use a clover- improved Wilson formulation |23| . In the following we will 



denote our simulation with the tree level clover coefficient csw = 1 as Run 
A and the simulation with tadpole- improved clover fermions, csw = l/^o^ 
with Run B. For k values and other details of the light quarks see Table |I[ 
The light quarks in Run A are rotated [Z3]: 

q(x) (l-~7-D^q(x), 

q(x) -> q(x) (l + -J- d) . (21) 



For the light quarks in Run B we use the normalisation ^Jl — 6k. 

The clover improvement removes lattice spacing errors at 0(a), and we 
expect the remaining leading errors for light quarks at zero momentum to be 
0(a s aAgcD) and 0(a 2 K 2 Q CD ). If we use for Aq C d a value around 300 MeV 
(aA v = 0.169 for our configurations), and a^(l/a) = 0.247, we estimate 
these errors to be ~4% and ~3% respectively. 

We use quenched gauge configurations at /3 = 6.0 on 16 3 x 48 lattices, 
generated by the UKQCD Collaboration. The configurations were fixed to 
Coulomb gauge. In Table |l| we list the ensemble sizes and the lattice spacings 
from light spectroscopy. Degenerate pion and rho masses, lattice spacings 
and results for K cr it and k s are taken for Run A from Ref. H, and for Run B 



from Ref. [2l|. The heavy quark parameters for both runs are given in the 
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same table. The variation of the ensemble sizes in Run B for different k and 
ttlq values between 45 and 67 is due to limited computer time. One has to 
note that for quenched configurations ratios of physical quantities generally 
differ from the corresponding ratios in the real world. Thus the values one 
obtains for the lattice spacing from different physical quantities are in general 
different. Averaging results from the IP — IS and 2S — IS splitting of the 
T, one obtains a -1 = 2.4(1) GeV at (3 — 6.0 Probably for heavy-light 
systems the appropriate lattice spacing is closer to the one determined from 
light hadron spectroscopy, and in the following we will use a~ l = 2.0(2) GeV 
to convert lattice results into physical units. This encompasses a -1 from m p 
for both runs (see Table [l]) and also lattice spacings from (see e.g. p6|, |27||) 
as well as from gluonic quantities |28 at (3 = 6.0. The heavy quark masses 



am'Q = 1.71, 2.0, 2.5 are in the region of the b quark. In our later simulation 

(Run B) we choose one heavy mass as high as arrLq = 8.0, because in Run 
A with the heaviest mass being 4.0 we find it difficult to extrapolate to 
the static limit. The value arrLQ = 1.71 corresponds to the b quark mass 
as determined from T spectroscopy with a~ x = 2.4(1) GeV |55|| . If one 



chooses the lattice spacing from light spectroscopy, a 1 ~ 2 GeV, arnq = 2.0 

gives approximately the same dimensionful b quark mass as the value arriQ 
= 1.71 used in the T simulations at a -1 = 2.4 GeV. This ambiguity in 
fixing lattice bare quark masses is typical of problems with the quenched 
approximation. 

At the source we use the following interpolating fields for the mesons: 

£ Q^x^ix.-xM^), (22) 



where T(x\ — X2) factorizes into a smearing function <f)(r = \x\ — X2I) and one 
of the 2x4 matrices in spinor space shown in table We calculate heavy- 
light current matrix elements with the smearing function <p being either a 
delta function or a ground or excited state hydrogen-like wave function [3U 



at the source and a delta function at the sink, and meson correlators with 
all combinations of smearing functions at source and sink. For the ground 



state, we use 



exp(-r/r ), (23) 
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Run A 


Run B 




Light quarks 


csw 




1.0 






1/rig 




rotated 




yes 






no 




K values 


0.1432 


0.1440 


0.1370 




0.1381 


am w 


0.386( + _ 


t) 


0.311( + _ J) 


0.4137( t 


l l) 


0.2940( 


am p 


0.51( t 


I) 


0.47( 1 ^ 


0.538( + 


I) 


0.463( + J) 


^crit 


0.14556(6) 




0.13926(1) 






0.1437( 1 


^ t)(m K ) 


0.13758(5) (m K ), 0. 13726(9) (m ) 


a~ x {m p ) 


2.0(2) GeV [ 


1 


1.99(4) GeV [ 


0] 




Heavy quarks 


(o) 
am> Q ' 


1.71 


2.0 


2.5 


4.0 


1.71 


2.0 


4.0 


8.0 


n 


2 


2 


2 


2 


2 


2 


1 


1 




Configurations 


number 


35 + time reversed 




45 - 


- 67 





Table 1: Simulation details. 




Table 2: Spin operators for mesonic states. 
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and for the radially excited state, 

4>{r) = Q) 3/2 (l - ^) exp(-r/(2r )), (24) 
choosing ar = 3. The smearing is applied on the heavy quark. 

3 Fitting procedure 

In this section, the lattice spacing a is set to 1. On the lattice, the decay 
constant can be extracted from the matrix element of the local current J^: 

fM={0\J k \B), (25) 

with k = 5 for pseudoscalar and k = 1,2,3 for vector decay constants. M 
denotes the meson mass. The correlation function of this current decays for 
sufficiently large times exponentially: 

C LL (t) - Z 2 L e- E ^ (26) 

We are using the notation C rs for the correlation functions, where the index 
r denotes the smearing function at the source, and the index s the smearing 
function at the sink. L stands for a delta function, 1 for the ground state 
and 2 for the excited state smearing function. E s[m is the bare ground state 
binding energy and Z L is related to the matrix element as follows: 

Zi = 7H? <01 JklBh (27) 

so that 

fy/M = V2Z L . (28) 

To distinguish between the uncorrected current and the current containing 
the first order correction in the i/m^-expansion, we will use the definitions 

(fVMr ncorr = ^=(0\4 0) \B), (29) 
v M 



and 

fVM^-^iOKJ^ + J^lB). (30) 
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The current corrections we denote as 



5{fVM)^^={Q\{J^\B). (31) 



Local meson operators, however, overlap considerably with excited states 
with the same quantum numbers. We therefore smear the heavy quark field 
at the source. The amplitude of the corresponding correlation function then 
contains also the density Z$, which we determine separately from smeared- 
smeared correlators. For sufficiently large times one has: 

C xl - A u e- E ^\ A n = Zl (32) 
C 1L -> A lL e- E ^\ An = Z S Z L . (33) 

We extract the matrix element Zl by fitting C\l and C\\ simultaneously to 
a single exponential, using a bootstrap procedure as described in |3l[] . Zl is 
calculated as the correlated ratio Ail/ y/ An. 

3.1 Run A 

For t m i n < 5, no single exponential fit to both C\l and C\\ is possible. 
After t m i n = 5 the amplitudes An and Ail and the ground state energy still 
decrease slightly, until t min is moved out to 9. The dependence of results 
on the fit interval is shown for m^' = 4.0 in Table [3[ For k = 0.1440 the 
decrease may be by more than la. We use a fit interval t m i n /t max = 9/25 
for all k and heavy mass values. Given that the statistics are poor, we have 
to discard the smallest eigenvalues of the covariance matrix in the singular 
value decomposition (SVD) algorithm. The fit results for the amplitudes 
from Run A are shown in Table f|. With our limited statistics, our bootstrap 
procedure generates certain ensembles on which multi-exponential fits to two 
correlators fail, but with the original ensemble of correlation functions we can 
do a simultaneous fit of C\l and Cil to 2 exponentials. This gives a value 
for the ground state energy which tends to be a little higher, but in general 
still compatible within one standard deviation with the results fitted with 
just one exponential The plateau of the static correlation functions 



sets in around tmin = 3 and persists for about 10 timeslices, which indicates 
that our ground state smearing functions work well for the static case. 

We determine the ratio of the axial current correction to the uncorrected 
axial current from a fit of the bootstrap ratio of the correlators of (0| J5 \B) 
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k = 0.1432 


k = 0.1440 


train/ tmax 


E ■ 

-'-'sim 


An 


An 


E ■ 

- LJ smi 


An 


Ail 


5/25 


0.516(6) 


211(12) 


2.37(10) 


0.497(6) 


202(11) 


2.15(11) 


6/25 


0.515(5) 


210(11) 


2.37(9) 


0.495(7) 


197(13) 


2.11(13) 


7/25 


0.514(6) 


210(13) 


2.34(14) 


0.493(8) 


194(12) 


2.05(14) 


8/25 


0.514(7) 


209(15) 


2.30(17) 


0.491(7) 


191(13) 


2.00(15) 


9/25 


0.513(7) 


210(15) 


2.28(19) 


0.491(8) 


189(13) 


1.97(16) 


10/25 


0.511(5) 


209(14) 


2.26(18) 


0.489(7) 


185(14) 


1.91(17) 



Table 3: Dependence of fit results from Run A at itLq = 4.0 on the fit 
interval. All quantities are in lattice units. 





k = 0.1432 


k = 0.1440 




(0) 
m Q 


An 


Ail 


An 


Ail 




1.71 


207(8) 


2.13(11) 


197(8) 


1.91(11) 




2.0 


206(8) 


2.17(10) 


195(14) 


1.96(15) 




2.5 


212(7) 


2.33(10) 


193(15) 


2.03(17) 




4.0 


216(12) 


2.59(17) 


199(14) 


2.27(19) 




static 


222(6) 


3.49(11) 


219(10) 


3.29(14) 




1.71 


202(16) 


1.66(9) 


188(18) 


1.48(19) 




2.0 


208(13) 


1.81(10) 


189(15) 


1.58(11) 




2.5 


210(11) 


1.99(12) 


188(15) 


1.71(13) 




4.0 


210(15) 


2.28(19) 


189(13) 


1.97(16) 


t(0) 
J k 


1.71 


215(20) 


2.02(18) 


208(16) 


1.83(15) 




2.0 


216(7) 


2.10(9) 


201(9) 


1.87(10) 




2.5 


218(9) 


2.24(13) 


202(5) 


1.98(10) 




4.0 


226(24) 


2.60(10) 


209(8) 


2.29(15) 



Table 4: Fit results for An and An from Run A. The columns on the left 
indicate the operators included in the currents. 
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Figure 1: Effective amplitudes for smeared-local (left) and smeared-smeared 
correlators (right) from Run A; k — 0.1432. The upper row shows the 1 S 
case at itlq = 2.0 with the 1/mg' correction to the current included, and 
the lower row, the static case. We choose the same scale of the x axis for the 
static correlators, to demonstrate the difference in the length of the plateau 
between NRQCD and static heavy quarks. 
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and (0| J5 \B) to a constant. As fitting interval we choose t min /t max = 9/25. 
Results are given in Table §. 

3.2 Run B 

The correlators from Run B are slightly noisier, as shown in the examples in 
Fig. |^. Also here we discard the lowest eigenvalues of the covariance matrix 
in our SVD inversions. For k = 0.1370, we choose t m i n = 6, for k = 0.1381, 
tmin = 5. The upper end of the fitting window is, for = 1.71, 2.0, and 

4.0, set to t max = 20, since for larger times the signal disappears. For toq 1 
= 8.0, the signal disappears earlier, so we choose for k = 0.1370, t max = 18 
and for k = 0.1381, t max = 15. The goodness-of-fit value Q of the fits is 
generally low (Q < 0.1). The fit results for the amplitudes with and without 
the current corrections are shown in Table |5|. 

At rriQ = 1.71, 2.0, and 4.0, the Cn correlation functions from Run B 
have after t ~ 14 (« = 0.1370) or t ~ 12 (k = 0.1381) a wiggle of ~ la. If 
we choose t max for the fits not to include this wiggle, we obtain a higher Q 
(~ 0.3 — 0.4), but the results for E sim and A\l are up to ~2<r smaller than the 
values listed in Table |^. We choose to extract the decay constant using the 
larger fitting range. However, this variation of the result with the fit range 
indicates that there is a fitting uncertainty in E s - im and A\l, associated with 
the choice of the fitting range, which could be for Wq' = 1.71, 2.0 and 4.0 
about twice as large as the bootstrap errors. This propagates into a fitting 
uncertainty of < 2a for the pseudoscalar matrix elements, and of ~ 2 — 3a 
for the vector matrix element (for chirally extrapolated light quarks it is in 
both cases < la). In the tables and figures, we give pure bootstrap errors 
on the results from Run B. We will include the fitting uncertainty where we 



quote our final results for decay constants (subsection |4.3|) . Note that the 
fitting uncertainty for the static case and for tuq = 8.0 cannot be estimated 
in a similar way since the signal from the correlator disappears much earlier. 

In a similar way to Run A, the plateau in the static correlation functions 
from Run B sets in slightly earlier and is clearly shorter than for the NRQCD 
correlators. For k = 0.1370 we use the fit interval t min /t max = 3/10 and for 
k = 0.1381, t m i n /t max = 3/9. Effective amplitude plots are shown in Fig. §. 

We also determine the ratio of the current corrections to the uncorrected 
current from the bootstrap ratio of their smeared-local correlators. For the 
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Figure 2: NRQCD correlators from Run B at m^' = 2.0 and k = 0.1370 with 
smeared source and local sink (left), and smeared source and sink (right). The 
upper row shows effective amplitudes, the lower row, effective masses. The 
smeared-local correlators include the 1/mg' correction to the current. 
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K 


= 0.1370 




= 0.1381 




(0) 
m Q 


E ■ 

-'-'sim 


sin 


rt-lL 


E ■ 

-'-'sim 


4-,., 


4-, r 


t(0) 


1 ti 


r\ r A ( A \ 

0.524(4) 


zUo( / ) 


z.Uo(/ ) 


0.503(5) 


nnn/ *7\ 
Z[)y{ 1 ) 


i.yy^y j 




2.0 


U.ozy^4 ) 


205(6) 


2.12(7) 


U.oUo^o ) 


208(8) 


2.02(9) 




4.0 


0.540(5) 


203(7) 


2.32(10) 


0.516(6) 


213(5) 


2.25(8) 




8.0 


0.540(6) 


204(9) 


2.54(13) 


0.511(8) 


204(7) 


2.39(11) 




static 


0.540(6) 


220(5) 


2.98(8) 


0.517(7) 


213(9) 


2.76(9) 


7(0) , r(l) 
J 5 -t- J 5 


1.71 


0.525(4) 


212(7) 


1.84(6) 


0.503(5) 


210(6) 


1.72(7) 




2.0 


0.529(4) 


214(7) 


1.94(6) 


0.503(6) 


208(8) 


1.76(7) 




4.0 


0.541(5) 


218(5) 


2.23(8) 


0.516(6) 


213(6) 


2.11(8) 




8.0 


0.542(8) 


218(8) 


2.62(8) 


0.513(8) 


211(8) 


2.37(11) 


t(0) 


1.71 


0.546(5) 


218(8) 


2.02(8) 


0.526(7) 


215(3) 


1.95(9) 




2.0 


0.548(5) 


218(6) 


2.10(7) 


0.519(7) 


214(5) 


1.93(8) 




4.0 


0.552(6) 


217(7) 


2.35(10) 


0.524(7) 


215(5) 


2.19(8) 




8.0 


0.550(6) 


219(6) 


2.62(10) 


0.518(10) 


216(3) 


2.49(9) 


4°' + 4" 


1.71 


0.544(4) 


220(5) 


2.14(6) 


0.522(6) 


216(7) 


2.00(8) 




2.0 


0.546(4) 


220(3) 


2.21(5) 


0.519(6) 


214(7) 


2.02(8) 




4.0 


0.550(5) 


220(6) 


2.46(9) 


0.528(7) 


216(4) 


2.29(8) 




8.0 


0.547(6) 


221(5) 


2.67(8) 


0.518(8) 


216(3) 


2.51(9) 



Table 5: Fit results for binding energies, A n and A 1L from Run B. The 
column on the left indicates the operators contributing to the currents. All 
quantities are in lattice units. 
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ratios at k = 0.1370 we choose the same fit intervals as described in the 
previous paragraph for the correlation functions. The results are shown in 
Table 0. The ratios at k — 0.1381 plateau later than the correlation functions, 



thus we use t min = 6 instead of 5. 
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Figure 3: Static effective amplitudes from Run B at k = 0.1370. 



4 Results 



We set a = 1 in this section, except in subsection 4.3 



4.1 Binding energies and meson masses 

In NRQCD, the heavy-light meson mass differs from the exponential falloff 
of the meson correlators, E S [ m , by the mass shift: 

A HL = M - E sim , (34) 

which depends on the renormalised heavy quark mass and the zero point of 
the non-relativistic heavy quark energy, Eq\ 

Ahl = Z m m>Q ~ Eq. (35) 

Z m is the heavy quark mass renormalisation constant. The meson mass can 
be determined nonperturbatively from the ratio of finite momentum and zero 
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Run A 


Run B 




k = 0.1432 


k = 0.1440 


k = 0.1370 


k = 0.1381 


(0) 






1.71 


0.1964(6) 


0.1986(5) 


0.1372(7) 


0.1438(7) 


2.0 


0.1696(2) 


0.1713(4) 


0.1215(5) 


0.1239(6) 


2.5 


0.1379(3) 


0.1389(4) 






4.0 


0.08800(13) 


0.08863(4) 


0.06561(18) 


0.0672(3) 


8.0 






0.003368(8) 


0.03447(11) 


(0) 

m Q 


\8(fVM) v /(fVM)^ ncorr \ 




1.71 






0.0455(2) 


0.0468(3) 


2.0 






0.04032(17) 


0.0410(3) 


4.0 






0.02178(5) 


0.0223(7) 


8.0 






0.01118(3) 


0.01151(3) 



Table 6: Fit results for the ratio of the current correction to the uncorrected 
matrix element, in lattice units. 5{f\fM)v has not been calculated directly 
in Run A. 
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momentum correlation functions [|32|. The dispersion relation of the heavy 
meson is given by: 



E sim (P) - ^sim(O) = yp 2 + M 2 — M. (36) 
We use the non-relativistic expansion of this: 

E S UP) = E sim (0) + {j^, (37) 

where M is the meson mass we want to determine. In this study the errors on 
the finite momentum correlators are rather large and we resort to different 
methods to calculate the meson mass. The shift can also be determined 
perturbatively [Q, or nonperturbatively from the mass shift (defined 
analogous to Ahl) in heavy- heavy systems: 

A HH = 2A HL (38) 

In Run B, the mass shifts for rriq = 1.71 and 2.0 are taken from the heavy- 
heavy results in Ref. [19], the shift for = 4.0 is from a heavy-heavy 



simulation which was part of this project. For rrLq = 8.0 we use perturbation 
theory since we expect the discretisation errors in heavy-heavy NRQCD for 
such large ttlq to be large. The determination of the binding energies and 



the meson masses from Run A is discussed in Ref. [31 



As shown in Table [7|, at K crit the pseudoscalar binding energies from Run 
A are ~ 10% lower than from Run B, the vector E sim are ~ 5% lower and 
the spin averaged E sim , ~6 — 8% lower. At « 8 , all the binding energies are 
~4 — 5% lower than the corresponding values from Run B. This difference 
corresponds to ~4cr. In the chirally extrapolated case error bars are larger 
and the difference between the two runs amounts to 2 — 3a. We expect Ahl 



to entirely depend on the heavy quark parameters [31], thus the difference in 
E sim translates directly into a difference in the meson mass. The values for 
Ahl used in Run B and the corresponding meson masses with chirally ex- 
trapolated light quarks are given in Table |8]. We note that the pseudoscalar 
and spin averaged masses from Run B seem to be enhanced by at least as 
much as the vector meson masses, which means that the enhancement of the 
clover term by tadpole improvement does not increase the hyperfine splitting 
by an amount greater than the statistical error. We came to the same con- 
clusion when we extracted the hyperfine splitting from the ratio of the 3 Si 



18 





^crit 


K s (m K ) 




Run A 


Run B 


Run A 


Run B 


(0) 

m n 

Q 




1.71 


0.437(19) 


0.480(12) 


0.489(6) 


0.513(3) 


2.0 


0.437(11) 


0.475(13) 


0.491(6) 


0.515(4) 


2.5 


0.439(12) 




0.496(7) 




4.0 


0.446(15) 


0.489(16) 


0.499(6) 


0.528(4) 


8.0 




0.482(16) 




0.527(6) 


oo 


0.461(15) 




0.511(6) 




(0) 
m Q 




E ■ ( 

- L -'sim V 


V) 




1.71 


0.464(15) 


0.498(15) 


0.509(6) 


0.532(4) 


2.0 


0.463(12) 


0.489(13) 


0.510(6) 


0.532(4) 


2.5 


0.464(9) 




0.512(6) 




4.0 


0.478(12) 


0.505(15) 


0.515(6) 


0.539(4) 


8.0 




0.488(17) 




0.532(6) 


(0) 
m Q 


E ■ 


1.71 


0.457(14) 


0.493(10) 


0.504(6) 


0.528(3) 


2.0 


0.456(12) 


0.486(13) 


0.505(6) 


0.527(4) 


2.5 


0.457(9) 




0.508(6) 




4.0 


0.462(13) 


0.501(13) 


0.511(6) 


0.536(3) 


8.0 




0.487(18) 




0.530(5) 



Table 7: Bare ground state energies from Run A and Run B at the strange 
and the chirally extrapolated light quark mass, in lattice units. 
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o) 

m Q 




M PS 


1.71 


1.73(10) 


2.21(11) 


2.0 


2.02(9) 


2.50(10) 


4.0 


4.07(5) 


4.56(7) 


8.0 


7.63(16) 


8.11(18) 



Table 8: Mass shifts used to calculate the meson masses in Run B, and 
chirally extrapolated meson masses using these shifts. The error on the per- 
turbative shift at itiq = 8.0 is an estimate of the contributions from higher 
orders in perturbation theory, obtained by squaring the one-loop contribu- 
tions to Eq and Z m . All numbers are in lattice units. 



and 1 S correlators (see Refs. |JT], |34|). Presumably the tadpole improvement 



of the light fermion affects the discretisation errors in the kinetic energy and 
the spin- magnetic energy of the heavy quark by the same fractional amount. 
Since the kinetic energy is much larger than the spin-magnetic energy, the 
spin-independent part of E sim receives the larger absolute shift due to tadpole 
improvement. 



4.2 Bare lattice decay matrix elements 

In this section we discuss the bare, unrenormalised lattice matrix elements. 



Renormalisation constants will be dealt with in subsection [O. 

The 1/M corrections in NRQCD can be separated into contributions of 
the kinetic and magnetic operator in the Lagrangian and the correction to 
the local current: 

fVM = (fVMr + — + + dM ° corr/6 ) . (39) 

J \ M M Ml 



The notation is chosen to be consistent with Refs. |3q, [I8| . For the axial 
current, du = 3 and for the vector current, du = — 1- With (fVM) 00 we 
denote the static matrix element. 
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4.2.1 fVM from Run A 



Decay matrix elements in lattice units from Run A, for the pseudoscalar, 
vector and spin averaged cases, are shown in Table || To obtain an estimate of 
the physical matrix elements and to be able to compare with other methods, 
we chirally extrapolate in the light quark mass to K cr i t and interpolate to the 
strange light quark mass at k s . These extrapolations should however be used 
with some caution, because we only have two k values for the light fermions. 
For comparison, we also list results which do not include the corrections 
to the current, denoted as {f \fM) uncorr . The current corrections from Run 
A are given in Table 0. On the same set of configurations, the matrix 
elements have been calculated using clover heavy quarks (csw = 1) by the 
UKQCD Collaboration j|, [36|]. Thus we are able to make a direct comparison 
between NRQCD and clover heavy quarks. Fig. |] shows (f\ /r M)ps as a 
function of the inverse pseudoscalar meson mass Mp$, both for clover and 
NRQCD (with the current correction included) heavy quarks. Note that 
for the mesons with clover heavy quarks we use M2, the 'dynamical' meson 
mass determined from the dispersion relation of the heavy meson [ fj6|j . The 
NRQCD results from Run A can be fit to a linear function in the inverse mass; 
a possible quadratic dependence of these results on the heavy mass cannot be 
resolved. The bare lattice matrix elements from both types of heavy quarks 
agree in the B region within errors. This is not necessarily true for the 
renormalised matrix elements. Renormalisation decreases {f\fM)ps from 
NRQCD in the B region by ~ 10%, as detailed in subsection £l|, whereas 
the current renormalisation constants have not been calculated for clover 
quarks with large mass. 

Neither of the results presented in Fig. [| extrapolates in the infinite mass 
limit to the static simulation result. Our static point agrees with UKQCD 
using different (Jacobi) smearing ||, and it appears that there is no sizeable 
excited state contamination which could artificially enhance the static result 
(see Table [U] and Fig. H). Instead, the reason seems to be that the quality 
of the results from Run A does not allow us to predict the slope at infinite 
quark mass correctly. In NRQCD calculations which include larger masses 
and better ensembles the results do extrapolate to the static point. This 
is detailed in the discussion of Run B in this paper (see Fig. |7]) and in 
Refs. |T?| In these studies we find that in the region of the B there is 
a quadratic contribution to the slope of fyM. We do however obtain good 
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k = 0.1432 


k = 0.1440 


^crit 


K s (m K ) 


(0) 
" l Q 




(fy/M) 


ancorr 
PS 




1 71 

1. 1 I 


0.209(8) 


0.193(8) 


0.160(18) 


0.199(7) 


9 n 


0.214(8) 


0.198(9) 


0.167(19) 


0.205(8) 




0.226(7) 


0.206(10) 


0.17(2) 


0.214(9) 


4.0 


0.249(11) 


0.227(13) 


0.19(3) 


0.235(12) 


oo 


0.332(7) 


0.315(9) 


0.281(19) 


0.321(7) 


(0) 
" L Q 


(f\/M) PS 


1 71 

J. . I J. 


0.165(7) 


0.153(8) 


0.128(21) 


0.157(6) 


2.0 


0.178(7) 


0.162(8) 


0.132(17) 


0.168(7) 


2.5 


0.195(7) 


0.176(9) 


0.138(19) 


0.183(8) 


4.0 


0.222(13) 


0.202(11) 


0.163(26) 


0.209(11) 


(0) 
" b Q 




(fy/M) 


uncorr 
V 




1.71 


0.195(10) 


0.180(9) 


0.15(2) 


0.186(8) 


2.0 


0.202(7) 


0.186(7) 


0.155(16) 


0.192(7) 


2.5 


0.215(9) 


0.197(8) 


0.161(13) 


0.204(8) 


4.0 


0.244(8) 


0.224(11) 


0.185(17) 


0.233(10) 


(0) 


(fVM)v 


1.71 


0.223(9) 


0.205(9) 


0.171(20) 


0.212(8) 


2.0 


0.226(8) 


0.210(9) 


0.177(19) 


0.216(8) 


2.5 


0.237(8) 


0.216(11) 


0.176(24) 


0.224(10) 


4.0 


0.256(11) 


0.234(13) 


0.191(24) 


0.242(12) 


(0) 


fy/M 


1.71 


0.199(9) 


0.183(8) 


0.153(18) 


0.189(7) 


2.0 


0.205(6) 


0.189(7) 


0.158(14) 


0.195(6) 


2.5 


0.218(8) 


0.199(7) 


0.163(11) 


0.206(8) 


4.0 


0.245(8) 


0.225(11) 


0.185(18) 


0.233(10) 



Table 9: Decay matrix elements fy/M in lattice units from Run A. 
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Figure 4: Unrenormalised matrix elements (fyM)ps in lattice units with 
NRQCD heavy quarks from Run A (circles) and with tree-level clover heavy 
quarks (diamonds) @, g§] at k = 0.1440. The cross denotes the static point 
from Run A. The lines are correlated fits; the errors on the extrapolations 
to infinite mass are for clarity slightly shifted from the origin. The arrow 
denotes the position of the B s meson in this plot (the light quark mass at 
k = 0.1440 is close to the strange quark mass). 
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Figure 5: Dependence of the static matrix element from Run A, in lattice 
units, on the fit interval. On the left, k = 0.1440, on the right, k = n CTit . 




Figure 6: Unrenormalised decay matrix elements from Run A, in lattice 
units. On the left, pseudoscalar (bursts) and vector (pluses) matrix elements 
without the current correction included. On the right, pseudoscalar (circles), 
and vector (squares) decay constants with the current corrections. The dia- 
monds on the right denote the spin-averaged case. The lines are correlated 
fits including the static point (cross), k = 0.1440. 
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re = 0.1432 


re = 0.1440 


recrit 


K s (m K ) 


(0) 
m Q 




1.71 
2.0 
2.5 
4.0 


-0.0412(16) 
-0.0364(13) 
-0.0312(10) 
-0.0219(9) 


-0.0383(17) 
-0.0340(15) 
-0.0287(14) 
-0.0202(11) 


-0.033(4) 
-0.029(3) 
-0.024(4) 
-0.017(2) 


-0.0394(14) 
-0.0349(14) 
-0.0297(13) 
-0.0208(11) 



Table 10: Current corrections to the decay matrix elements from Run A, in 
lattice units. 



t min /tmax 


re = 0.1432 


re = 0.1440 


recrit 


re s 


3/10 


0.333(6) 


0.319(6) 


0.293(8) 


0.324(5) 


4/10 


0.334(7) 


0.317(7) 


0.285(9) 


0.323(7) 


5/10 


0.333(9) 


0.316(9) 


0.282(15) 


0.322(8) 


6/10 


0.337(10) 


0.319(13) 


0.29(4) 


0.326(11) 


7/10 


0.343(13) 


0.324(13) 


0.29(3) 


0.331(14) 



Table 11: The static fyM from Run A for various fit intervals, in lattice 
units. 

quadratic fits of the results from Run A if we include the static point into 
the fit. In Figure |6| we show the NRQCD and static decay matrix elements 
from Run A at re = 0.1440 with a correlated fit to a quadratic function in the 
inverse meson mass that includes the static point. For the results without 
inclusion of the current corrections, at the left in Fig. |5|, we find Q = 0.94 
for the pseudoscalar and Q = 0.96 for the vector matrix element. The vector 
matrix element is smaller than the pseudoscalar matrix element. On the 
right, we plot the matrix elements with the correction to the current. The 
current correction to (fyM)y has not been simulated in Run A, but we 
estimated the corrected vector current (see Table [|) using the axial current 
corrections in Table |] and the relation: 
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which follows from Eq. (^) with the appropriate values for du- We find 
(fv~M)ps in the B region to be much lower, around 50%, than in the static 
case. The current correction gives a contribution of relative size 22 % in the 
B region. The spin averaged matrix element, defined as: 

1 r 



fVM = - [(fVM)p S + 3(/VM) y ] , (41) 

is also shown in Fig. [| It can be calculated from the uncorrected matrix 
elements, since the current correction drops out after spin averaging (see 
Eq. fl3"9])). With the current corrections included, the quadratic fit of the 
pseudoscalar matrix element gives Q = 0.58, of the vector matrix element 
Q = 0.95, and the spin averaged, Q = 0.93. 

4.2.2 fVM from Run B 

In Fig. ^ we show the pseudoscalar, vector and spin averaged decay matrix 
element from Run B with correlated fits of the bare simulation results to a 
quadratic function in the inverse pseudoscalar meson mass. The results of 
these fits for the pseudoscalar matrix elements for both light k values are 
shown in Table |nj Within errors, the infinite mass limit of the NRQCD 
matrix elements are in agreement with the static simulation result. We find 
the relative slope of fyM to be of the order of 2.5 GeV. The error on the 



slope is however large (see Table 14), because of our small ensemble size, 



and the number of degrees of freedom for these fits is low. Our results 



from a simulation using dynamical configurations with higher statistics [21 
indicate that the slope actually gets smaller with decreasing light quark mass. 
Renormalisation decreases the NRQCD decay matrix elements relative to 
the bare ones (see subsection |4.3|) , but their infinite mass limit is still in 
agreement with the (renormalised) static matrix elements. 

To study the behaviour of the lattice matrix elements for the current 
corrections in the infinite mass limit, we perform correlated fits of 8(fy/M) ps 
and S(f\/M)v as functions in the inverse pseudoscalar meson mass. Fits to 
a linear function have Q -C 1, but fits to a quadratic function work well. The 



results are presented in Table [15|. We find that the infinite mass limit of our 
results is in a reasonable agreement with zero. In Fig. the axial and vector 
current corrections for k = 0.1370 are shown. 

From the heavy mass dependence of the spin average of the matrix ele- 
ments, we can, using Eq. (|39|), extract the contribution of the kinetic energy 
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0.0 0.2 0.4 

l/Mp S 



Figure 7: Unrenormalised decay matrix elements from Run B at k — 
0.1370 plotted against the inverse pseudoscalar meson mass. Circles de- 
note (fVM)pS: squares (fyJ~M)y, and diamonds the spin-averaged matrix 
element. Lines denote correlated fits; we indicate the error bar on the ex- 
trapolation of the pseudoscalar matrix element to infinite mass. The static 
point (cross) is slightly shifted from the origin for clarity. All quantities are 
in lattice units. 
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0.0 0.2 0.4 

l/Mp S 



Figure 8: Unrenormalised current corrections from Run B, in lattice units, 
at k — 0.1370, plotted against the inverse pseudoscalar meson mass. Circles 
denote —5{f\fM) PS , and squares 5(fy/M) v . The lines denote correlated fits 
of \5(fy/~M)\ to a quadratic function in 1/Mpg. Where not shown, error bars 
are smaller than the symbols. 
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k = 0.1370 


k = 0.1381 


^crit 


K s (m K ) 


(0) 
m Q 


(/v / M)™ r 


1 71 
1. / 1 


0.212(5) 


0.201(6) 




0.188(14) 


0.206(4) 


z.U 


0.219(5) 


0.203(4) 




0.187(12) 


0.211(3) 


a n 
4.U 


0.242(8) 


0.225(6) 




0.206(15) 


0.233(4) 


8.0 


0.262(9) 


0.243(8) 




0.223(18) 


0.252(6) 


oo 


0.293(6) 


0.274(7) 




0.254(15) 


0.283(5) 






(/V 


/ M) 


PS 




1 71 
1. / 1 


0.186(4) 


0.173(5) 




0.158(12) 


0.179(3) 


o n 


0.196(4) 


0.178(4) 




0.158(9) 


0.186(3) 


4.0 


0.232(6) 


0.210(6) 




0.186(14) 


0.221(4) 


8.0 


0.261(7) 


0.237(9) 




0.211(17) 


0.248(5) 


(0) 
m Q 


(/v / M)™ r 


1 71 


0.202(5) 


0.194(9) 




0.185(18) 


0.198(5) 


9 


0.210(5) 


0.192(6) 




0.174(12) 


0.201(4) 


4.0 


0.236(9) 


0.217(6) 




0.197(16) 


0.226(5) 


8.0 


0.261(6) 


0.240(9) 




0.218(19) 


0.250(6) 


(0) 




(A 


ru) y 




1.71 


0.213(5) 


0.198(6) 




0.181(16) 


0.205(4) 


2.0 


0.220(5) 


0.200(7) 




0.179(14) 


0.210(4) 


4.0 


0.244(8) 


0.226(7) 




0.206(16) 


0.235(5) 


8.0 


0.266(8) 


0.243(8) 




0.219(17) 


0.254(6) 


(0) 




/ 


s/M 




1.71 


0.207(5) 


0.191(5) 




0.175(12) 


0.199(4) 


2.0 


0.214(5) 


0.195(6) 




0.174(12) 


0.204(4) 


4.0 


0.241(7) 


0.222(7) 




0.201(13) 


0.231(4) 


8.0 


0.265(7) 


0.242(10 


) 


0.217(18) 


0.253(6) 



Table 12: Decay matrix elements from Run B, in lattice units. 
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k = 0.1370 


k = 0.1381 


^crit 


Kg(m K ) 


(0) 
m Q 


8(fVM)p S 


1.71 


-0.0291(7) 


-0.0287(9) 


-0.028(2) 


-0.0289(5) 


2.0 


-0.0265(6) 


-0.0252(6) 


-0.0238(16) 


-0.0258(4) 


4.0 


-0.0159(5) 


-0.0151(4) 


-0.0143(10) 


-0.0155(3) 


8.0 


-0.0088(3) 


-0.0084(3) 


-0.0079(6) 


-0.0086(2) 


(0) 

m o 


5(fVM)v 


1.71 


0.0092(2) 


0.0091(4) 


0.089(9) 


0.091(3) 


2.0 


0.0085(2) 


0.0079(2) 


0.073(5) 


0.0816(14) 


4.0 


0.00513(18) 


0.00485(13) 


0.045(4) 


0.0498(11) 


8.0 


0.00292(7) 


0.00277(10) 


0.026(2) 


0.0284(7) 



Table 13: Current corrections to the decay matrix elements from Run B, in 
lattice units. 



K 


ai 


a 2 /ai 


a 3 /ai 


0.1370 
0.1381 


0.30(2) 
0.28(3) 


-1.2(4) 
-1.4(7) 


1.6(1.5) 
2.6(2.4) 



Table 14: Results of correlated fits of (fy/~M) PS in lattice units from Run B 
to the function ai + a 2 /M PS + a 3 /Mp s . 







a 2 


a 3 


K 


5(fWM) PS 


0.1370 


-0.0016(13) 


0.091(11) 


-0.10(4) 


0.1381 


-0.0016(14) 


0.087(14) 


-0.09(4) 


K 


5(fVM) v 


0.1370 


-0.0016(13) 


0.091(11) 


-0.10(4) 


0.1381 


-0.0004(5) 


0.028(5) 


-0.033(15) 



Table 15: Result of correlated fits of the current corrections from Run B to 
the function a\ + a 2 /Mp S + a 3 /Mp s . All quantities are in lattice units. 
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Figure 9: Spin averaged decay matrix elements, unrenormalised, in lattice 
units, at k s , plotted against the inverse pseudoscalar meson mass. Squares 
are results from Run A, and circles, from Run B. Lines denote correlated fits 
to all the points including the static. 
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Run A 


Run B 


^crit 


-1.6(2) 
-2.0(5) 


-1.05(15) 
-1.2(5) 



Table 16: Results for Gkin from correlated fits of (JvM) to a quadratic 
function in 1/Mpg, in lattice units. 



to fVM pi 



fVM = (/VM)°° 1 + 



G 



kin 



M 



(42) 



We fit the spin-averaged matrix elements from Run A and Run B, including 
also the static point, to a quadratic function in 1/Mps- Results in lattice 
units for the strange and the chirally extrapolated light quark mass can be 
found in Table [0]. Fig. |9] shows f\fM at k s from both runs. For Run B, 
Gkin is found to be ~ 2a smaller than for Run A. However it has to be 
noted here that the fits include 5 points and thus there are only 2 degrees 
of freedom. For a definite conclusion it would be necessary to include more 
simulation results at heavy quark masses. Moreover, things are expected to 
change after inclusion of renormalisation constants. We expect the heavy 
mass dependence for each run to change after renormalisation in a different 
way. In particular, it appears that rotating the light quark in Run A with the 
7 ■ D operator introduces an additional, heavy mass dependent, contribution 
to Za and Zy in perturbation theory. 



4.2.3 

We can also study the behaviour of the ratios of axial and vector matrix 
elements in the heavy quark limit. The ratio of the matrix elements without 
the current correction should give an estimate for the contribution of the 

m 



spin-magnetic interaction in the Hamiltonian to j\JM in Eq. (|3 

r - (/VM) 



, uncorr 
* PS 



(fVM)f 



oc 1 + 



8G 



hyp 



M 



(43) 
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For the slope of the ratio of the matrix elements with the current correction 
one expects: 

R _ (fVM)pS (xl + % G hyp + 2G corr /3 ^ 



(fVM) v M 
The correlated ratios of axial and vector matrix elements from both runs 




0.7 



0.7 



0.0 



0.2 
l/Mp S 



0.4 



0.0 



0.2 
1/M pg 



0.4 



Figure 10: Ratio of unrenormalised axial and vector matrix elements at k = 
k s from Run A (left) and Run B (right). Bursts refer to results without the 
current correction included, squares to results with the current correction. 
The lines denote correlated fits to the ratios. 



are listed in table |17[ The ratios can for both runs be fit to a linear function 
in 1/Mps- The fit results are presented in Table [L^. As expected, we find 
the extrapolation to the infinite mass limit to be in good agreement with 
one. We also perform fits with the value of the ratio fixed to one at infinite 
mass. In some cases these fits are slightly worse (see Table [18); however 
the slopes from both fit methods agree with each other. We find the slope 
of the uncorrected ratio and thus Gh yp to be slightly larger for Run A than 
for Run B, but the difference is not statistically significant. We note that 
the uncorrected ratios are independent of the light quark mass, which we 
also found for the B* — B hyperfine splittings |^TJ. The absolute value of 
the slope of the ratio R turns also out to be larger in Run A than in Run 
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Run A 




« = 0.1432 


k = 0.1440 


^crit 


K s 


(0) 
m Q 


(/ 


v / M)^ cor 7(/v / M)r corr 


1.71 


1.07(5) 


1.07(5) 


1.07(18) 


1.07(4) 


2.0 


1.06(4) 


1.07(4) 


1.09(14) 


1.07(3) 


2.5 


1.05(3) 


1.05(3) 


1.05(20) 


1.05(3) 


4.0 


1.02(3) 


1.02(5) 


1.02(18) 


1.02(4) 


(0) 
m Q 


UVM) PS /{fVM) v 


1.71 


0.74(3) 


0.74(3) 


0.75(13) 


0.74(2) 


2.0 


0.78(2) 


0.77(3) 


0.76(14) 


0.78(2) 


2.5 


0.826(16) 


0.82(3) 


0.80(12) 


0.82(2) 


4.0 


0.87(4) 


0.86(4) 


0.86(19) 


0.86(4) 


Run B 




k = 0.1370 


k = 0.1381 


^crit 


K s 


(0) 
m Q 


(/ 


v / M)^ cor 7(/v / M)r corr 


1.71 


1.05(2) 


1.04(5) 


1.03(12) 


1.04(3) 


2.0 


1.041(16) 


1.057(19) 


1.08(5) 


1.049(12) 


4.0 


1.028(19) 


1.036(18) 


1.05(5) 


1.032(13) 


8.0 


1.002(18) 


1.010(17) 


1.02(5) 


1.006(13) 


(0) 


(fVM) PS /(fVM)v 


1.71 


0.87(2) 


0.87(4) 


0.88(9) 


0.87(2) 


2.0 


0.890(14) 


0.887(17) 


0.88(4) 


0.889(11) 


4.0 


0.951(14) 


0.93(3) 


0.90(8) 


0.94(2) 


8.0 


0.983(9) 


0.974(15) 


0.96(4) 


0.978(9) 



Table 17: Ratio of axial and vector matrix elements from Run A and Run 
B, in lattice units. 
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Run A 


Run B 


ai 


a 2 


Q 




a 2 


Q 



mncorr 



0.96(8) 
1 (fixed) 


0.25(19) 
0.16(5) 


0.97 
0.92 


0.998(18) 
1 (fixed) 


0.13(7) 
0.12(3) 


0.69 
0.86 


R 


1.0(5) 
1 (fixed) 


-0.55(13) 
-0.56(3) 


0.74 
0.90 


1.017(11) 
1 (fixed) 


-0.33(5) 
-0.27(2) 


0.93 
0.42 



Table 18: Results from correlated fits of the ratio of axial and vector decay 
matrix elements, in lattice units, to the function ai+G^/^fps 1 - k = k s . R uncorr 
is the ratio of the uncorrected matrix elements, R has the current corrections 
included. 



B. Since the sign of the combination 8Gh yp + 2G corr /3 is opposite to Gh yp , 
this means that the current correction in Run A is considerably larger than 
in Run B. Another, more precise, way to determine G corr is from the ratio 



of the current corrections to the uncorrected current [|T^]. From Eq fl39| ) it 
follows that 

G corr = hm p(mg>)ES lim H^W^W 

(o) rv V ' (o) ( t./M\ uncorr 

rrig'— >oo nig ^oo \Jv lvl )ps 

We fit pirriq) from Run A to a linear function in 1/Mps and find G corr = 
—0.739(5). If we include all values for p from Run B in a linear fit, we 
obtain a bad fit (Q = 0.01). The situation improves with a quadratic fit 
(Q = 0.16), but we obtain a better result from a linear fit, omitting the point 
at ttlq = 4.0 from the fit. This gives Q = 0.44, and for G corr we obtain 
—0.5703(16). Both fits are shown in Fig. O. In summary, we find each of 
the three different 1/M corrections to fvM, as well as the total (compare 
Tables || and 112, and Tables and |13|), in Run A to be larger than in Run 
B. 
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Run A 


Run B 


(0) 

rriQ 


K — K cr it 


K = K s 


K — K cr it 


K = K s 


1.71 


0.2030(16) 


0.1978(3) 


0.1494(18) 


0.1403(5) 


2.0 


0.1747(10) 


0.1707(3) 


0.1265(15) 


0.1227(3) 


2.5 


0.1410(8) 


0.1385(3) 






4.0 


0.08988(10) 


0.0884(3) 


0.0688(7) 


0.06642(17) 


8.0 






0.0353(3) 


0.03409(6) 



Table 19: Ratio \5(fVM)p S /(f\ r M)'f l s corr \ of the axial current correction to 
the uncorrected matrix element, in lattice units, at the chirally extrapolated 
and the strange quark mass. 




Figure 11: Ratios of unrenormalised current corrections at k — k s plotted 
as a function of the inverse heavy meson mass. Squares denote results from 
Run A, circles results from Run B; the lines are correlated fits to the function 
a x + a 2 /M PS . 
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4.2.4 f Bs /f Bd 

The ratio f B Jf Bd can be used for an extraction of the ratio of CKM matrix 
elements |14s/14d| [0- HQET predicts that it is up to corrections 0{m s — 
mdj/mQ independent of the heavy quark mass. The renormalisation con- 
stants are expected to be very weakly dependent on the light quark mass. 
We approximate the physical ratio fB s /fB d with the ratio of the unrenor- 
malised lattice matrix elements: 



fa ^ (NM) Bs 
f Bd (fyffl) Bd 

Note that the renormalisation constants used in the following subsection 
assume massless light quarks. The results for f B J ' f Bd are listed in Table EU 



As expected, they are within errors independent of the heavy quark mass. 
We find that the ratio is larger for Run A, although not significantly. This 
disagreement might be a reflection of the different discretisation effects in 
both runs. 



o 

m o 


Run A 


Run B 


1.71 


1.23(18) 


1.13(8) 


2.0 


1.26(15) 


1.17(7) 


2.5 


1.32(12) 




4.0 


1.25(25) 


1.18(8) 


8.0 




1.17(9) 


static 


1.13(6) 


1.10(6) 



Table 20: Ratios of decay constants f B J f Bd . 



4.3 Renormalised matrix elements 

The results presented in the previous subsection are at tree level, i.e. they 
do not include the renormalisation constants which are required to match 
between the matrix elements in the effective theory on the lattice and the 
matrix elements in full QCD in the continuum. 
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In NRQCD, the operators that contribute to the heavy-light current mix 
under renormalisation. The matrix elements of the bare NRQCD current 
operators on the lattice, j/*) t , the renormalised NRQCD current operators, 
J$?J n , and the heavy-light current in full QCD, Jqqd, are thus related by: 

(Jqcd) = E^(4:l) = £^£^V£>- (47) 



At 0{1/itlq), three operators contribute to the axial vector current: 



7(0) 

J latt 



qLQ, 



r(l) 
J latt 



1 



-' (oj9i2^ • DQ, J® 



2m 



2m 



(o) 



q{ 2 D -BQ. (48) 



(2^ (2) 

JiJ t does not contribute at tree level. Although we did not simulate J /ati , we 
know its matrix element exactly, since translation invariance implies that at 
zero momentum (J® t ) = 

Expanding the renormalisation constants rji and through 0(a), Eq. ([IT]) 
becomes 



(Aqcd) 



+ 



1 + a(B - Coo 
1 + a(B 1 - Coi 



+ Ks 2 -Co2-Ci 2 )](^). 



U-\{c Q + c q )) 

Cxx-\{C Q + C q )) 



( 7 (0) ) 

\ d lattl 

( 7 (1) ) 

\ J lattI 



(49) 



The coefficients Bi originate from the vertex and wave function renormali- 
sation in the continuum, the Qj denote the vertex renormalisations on the 
lattice, and Cq and C q the heavy and the light wave function renormalisa- 
tion on the lattice, respectively, a is calculated using the two-loop formula 
We assume a reasonable choice for the scale for ay lies be- 
1/a and q* = ir/a. The contribution of the rotation to the 



for ay pB | 
tween q* = 

heavy-light vertex with NRQCD heavy quarks has not been calculated, so 
we cannot include the renormalisation constants to Run A. The renormalised 
matrix elements for Run B have been calculated using Bi, Qj and Cq from 
Ref. p7| , and C q from Ref. []3~9[| . The results for chirally extrapolated and 



strange light quark masses are given in Table [21]. Our renormalised matrix 
elements include an 0(aa) discretisation correction to (jjatt)- The origin of 
this correction and its relation to Eq. ( |49"1) will be discussed in future publi- 
cations [57 , 4T| [4"2]j . In Fig. [12] we compare the renormalised a 3 ^ 2 (f\/M)ps 
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Figure 12: Matrix elements a 3 ^ 2 (f^/M)p S from Run B at n s . The circles 
denote unrenormalised matrix elements; the squares denote renormalised re- 
sults using aq* = 1, and the diamonds, using aq* = n. 
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(0) 

anriQ> 


^crit 




s 


aq* = 1 


aq* = Ti 


aq* = 1 


aq* = ti 


1.71 


0.143(14) 


0.149(14) 


0.160(4) 


0.166(4) 


2.0 


0.141(11) 


0.149(12) 


0.161(3) 


0.170(3) 


4.0 


0.146(13) 


0.163(14) 


0.167(3) 


0.186(4) 


8.0 


0.146(15) 


0.172(16) 


0.167(5) 


0.195(5) 


oo 


0.146(11) 


0.186(12) 


0.163(4) 


0.208(4) 



Table 21: Renormalised a?/ 2 {fyM)ps from Run B at strange and chirally 
extrapolated light quark masses. 

for both values of q* with the unrenormalised matrix element. After renor- 
malisation, the slope for large masses is remarkably smaller than before. For 
aq* = 1, it is within errors in agreement with zero. For aq* = tt, the M 
dependence of the matrix elements is approximately linear, and we estimate 
the relative slope (the slope divided by the value at infinite mass) to be of 
the order of —1 GeV. The latter value is roughly in agreement with previous 
calculations; e.g. in Refs. [§, |[ [7| the relative slope is ~— 1 GeV. Since the 
results at heavy masses vary considerably depending on q*, we do not want 
to make a more quantitative statement here. 

For fs, we use the result at arnq = 2.0, the point in our simulation 
whose mass is closest to the actual B meson (see Table ||): 

q* = l/ a : f B = 0.174(28)(26)(16) GeV 

/ B , = 0.198(8)(30)(17) GeV 
q* = Tx/a: f B = 0.183(32) (28) (16) GeV 

f Bs = 0.209(8) (32) (17) GeV (50) 

This is to be compared with the tree level result, f B = 0.195(22) (29) (7) GeV, 
and f Bs = 0.229(8) (35) (8). The first error bar is the statistical error, inflated 



by a factor of 2 to take the fitting uncertainty of ~ la (see Sec. |3.2|) into 
account. For f B , the statistical error gets in addition enlarged due to the 
chiral extrapolation. The second error bar stems from the uncertainty in 
the determination of a. The third one consists of the estimated error of the 
perturbative calculation and due to neglected orders in the 1/M expansion, 
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and was determined as follows: The uncertainty from the choice of q* is 
5%, which can be used as an estimate of the uncertainty from higher orders 
in perturbation theory. Another estimate of this uncertainty, 62, can be ob- 
tained from the relation a 2 ^ = (a5i) 2 , where Si is the one-loop contribution. 
Renormalisation decreases the decay constant for q* between ir/a and 1/a 
by 9 — 10% from the bare value. Squaring this yields a contribution of ~ 1% 
from higher orders. The numerical errors on the integrals in the perturba- 
tive calculation are estimated to be ~ 2 — 3% and propagate to the final 
result after being multiplied by a. Another source of error are higher order 
contributions in the 1/tuq expansion. A calculation that includes also the 

0{1/{itlq ) 2 ) corrections [41], HI indicates that these corrections are in the 



region of the B of the order of ~3 — 4%. Adding this to the errors from the 
perturbative calculation, we quote a systematic error of 9% of the average 
from both q* on our results in Eq. fl57}|), which is represented by the third er- 



ror bar. For the unrenormalised matrix elements, the third error bar consists 
only of the 0(1/ '(toq) 2 ) correction. The largest errors on our results on fg 
come therefore from the statistical and fitting error, which is magnified due 
to chiral extrapolation, and the uncertainty in a. Note that at the B, the 
error from the uncertainty in q* is much smaller than for the higher masses. 
The quenching error is partly reflected in the error we quote for the lattice 
spacing. There are indications that unquenching leads to a larger value for 



f B [0, 21 



For the static matrix element from Run B we obtain: 

q *= 1/a: /jf = 0.180(14)(27)(50) GeV 

q * = 7T /a: = 0.229(15) (34) (50) GeV, (51) 



while the bare result is a s/2 (fVM)°° = 0.254(15). An estimate of aq* = 2.18 
has been given in Ref. [f43j for static heavy and Wilson light quarks. In 



the static case the one-loop contribution to Za is larger than for NRQCD 
around the b quark mass, thus also the variation with q* is larger (24%). Here, 
the dominant error appears to originate from higher orders in perturbation 
theory. 

For the static matrix element, the renormalisation constant has been cal- 



culated also for rotated clover fermions j|4|, |37]; however the O(aa) dis- 



cretisation correction has not been determined for this case. Without this 
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discretisation correction, we obtain for the static fs from Run A: 

q* = l/a: a 3 / 2 (f^M) 00 = 0.154(10) /f = 0.190(12)(21) GeV 

q* = 7T /a: a 3 / 2 (f^M) 00 = 0.202(14) /~ = 0.249(17) (27) GeV, (52) 

compared to the bare result a 3 / 2 (/v / M)°° = 0.281(19). Note also that the 
result in Eq. (|52]) is not tadpole- improved. 

5 Summary and Conclusions 

We report on a study of quenched heavy-light decay constants with non- 
relativistic heavy quarks in a mass range around the b quark and heavier. 
Both the NRQCD Lagrangian and heavy-light current are correct through 
0(l/mg' 1 ). We performed two simulations, one uses clover light quarks with 
tree-level clover coefficient (Run A), the other uses tadpole- improved light 
quarks (Run B). 

We investigated the 1/M behaviour of the unrenormalised decay matrix 
elements fVM, M being the heavy-light meson mass. In the mass region 
of the B meson, the correction to the static limit is large (for Run A ~ 
50%, for Run B ~35 — 40%), before renormalisation constants are included. 
We disentangle the various 1/M corrections to the decay constants, and 
compare their size between Run A and Run B. The differences are small for 
the contributions of the hyperfme and kinetic term in the action, but sizeable 
for the current correction matrix element, at least at tree level. At the B, 
we find for the bare axial matrix element a difference of 18% (1 — 2a) for 
chirally extrapolated light quarks and 10% (2 — 3a) for strange light quarks. 
We expect this difference to be partly caused by a reduction of the 0(aa) 
errors in Run B due to tadpole improvement. The fact that we use a rotation 
with a derivative operator in Run A and a normalisation with \/l — 6k in 
Run B also introduces a difference in the discretisation effects between the 
two runs. 

We compare fy/M using NRQCD with results using clover (csw = 1) 
heavy quarks generated by UKQCD. Renormalisation constants are not in- 
cluded. In the region of the b quark both methods agree within errors. 
However, they behave quite differently at large masses, such that the clover 
results cannot be made to extrapolate to the static limit. 
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We calculated the renormalisation constant Za for NRQCD in one-loop 
perturbation theory, taking into account the mixing between the current 
operators. For the first time, we present renormalised pseudoscalar decay 
constants, and a value for /# from NRQCD where all the matching factors 
through 0(a/M) are included. The bare matrix elements show a larger slope 
in 1/M than the results from calculations with relativistic heavy quarks, but 
the heavy mass dependence of the renormalised matrix elements f\/~M is 
much milder than before renormalisation. 
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